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ABSTRACT 

We examine the properties of a morphologically-selected sample of 0.4 < z < 1.0 spheroidal galaxies 
in the GOODS fields in order to ascertain whether their increase in abundance with time arises 
primarily from major mergers. To address this question we determine scaling relations between 
the dynamical mass Mdym determined from stellar velocity dispersions, and the stellar mass M^, 
determined from optical and infrared photometry. We exploit these relations across the larger sample 
for which we have stellar masses in order to construct the first statistically robust estimate of the 
evolving dynamical mass fu nction (MF) oyer < z < 1. The trends observed match those seen in 
the stellar mass functions of iBundv et al.l (|2005f) regarding the top-down growth in the abundance of 
spheroidal galaxies. By referencing our dynamical masses to the halo virial mass, Myir, we compare 
the growth rate in the abundance of spheroidals to that predicted by the assembly of dark matter 
halos. Our comparisons demonstrate that major mergers do not fully account for the appearance of 
new spheroidals since z ^ 1 and that additional mechanisms, such as morphological transformations, 
are required to drive the observed evolution. 

Subject headings: galaxies: evolution — galaxies: elliptical and lenticular — galaxies: kinematics and 
dynamics — galaxies: mass function 



1. INTRODUCTION 

Substantial progress has been made in characterizing 
the primarily old, red, and passively-evolving stellar pop- 
ulations of spheroidal galaxies (defined here to include el- 
liptical and le nticulars) from z~l to the present day (see 
iRenzinil I2006L for a review). Their evolving abundance 
contains valuable clues to the processes by which they 
form, as well as insight into the history of cosmic mass as- 
sembly. In recent years, many surveys have converged on 
a broad evolutionary pattern in which it is claimed that 
galaxies with the largest stellar masses evolve predomi- 
nantly into spheroidals by z > 1 with little subsequent 
growth. Meanwhile, intermediate and lo wer mass sys- 
tems continue to form at later epochs fe. g..lBundy et ail 
20051120061: i Tanaka et all 120051: IF rancesc him et al.ll2006t 
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1006; Abraham et al. 2007), 

Some authors (e.g., Ivan DokkumI 120051) have ques- 
tioned the absence of growth among massive galax- 
ies, pointing to a high rate of tidal features in lo- 
cal examples , indicative of rec ent merging. It is like- 
wise argued (IDrorv et al.l 120041: Ivan der Wei et al] 120061: 



iMaraston et al.ll2006l : iKannappan k Gawiseill2007h that 
stellar masses derived from optical-infrared photometry 
may suffer biases that lead to large uncertainties in the 
inferred growth rate. Beyond such concerns, interpret- 
ing the observed evolution in the context of the ACDM 
paradigm necessitates appealing to potentially uncertain 
semi-analytic models. Clearly what is needed is a way 
to connect observations of galaxy growth directly to the 
assembly history of the halos in which they reside. 
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Spheroidal galaxies offer particular advantages in this 
regard. First, because of their regular mass profiles, spec- 
troscopic observations can be used to infer the dynamical 
mass, bypassing luminous quantities. Second, because of 
their pressure-supported dynamical configurations, they 
are e xpected to result, for the most part , from mergers 
(e.g.. lToomr3ll977HSpringel et al.ll2005af l. Accordingly, 
their evolving abundances can be compared directly to 
the expected assembly history of dark matter halos. In 
this Letter we develop a technique to estimate dynami- 
cal and halo masses for a large sample of field spheroidals 
and use it to interpret their evolution in the context of 
halo assembly as predicted in numerical simulations. We 
assume V,m = 0.3, flA = 0.7, Hq = 70/i7o km s~^ Mpc~^. 

2. OBSERVATIONS AND MASS ESTIMATES 

The dynamical mass (hereafter Mdyn) measurements 
us ed in this L e tter a. re based on the analysis presented 
in iTreu et al] ()2005D . to which we refer for details. 
Briefly, this sample comprises 165 field visually-classified 
spheroidals selected from ACS imaging in the northern 
GOODS field to a magnitude limit of zab < 22.5 and 
studied with the DEIMOS spectrograph on Keck II. A 
subset of 125 galaxies have reliable stellar velocity disper- 
sions, stellar masses, and surface photometry fits, which 
provi de B-hand effective radii, Res (see iTreu et al] 
l2005f) . Dynamical masses were derived using the for- 
mula, Mdyn = Kya^Res/G where is the virial coeffi- 
cient and a is the central velocity dispersion, obtained 
by increasing the observed dispersion by 10% to cor- 
rect to the standard circular aperture of diameter i?cB/4. 
Detailed modelling of fiber as well as 3D spectroscopy 
of nearby ellipticals supports the use of this equation 
(Padmanabhan et al. 2004; Cappellari et al. 2006) with 
little uncertainties due to the overall flatness of the ve- 
locity dispersion profile. 

With the kinematic sample defined in this way, the 
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investigation that follows is based on the much larger 
sample of galaxies with stellar masses (h ereafter M*) 
in bo th GOODS fields wfth zab < 22.5 (|Bundv et al.l 
l2005f) . We refer to that paper for further details but sum- 
marize the key features. For 0.4 < z < 1.0, spectroscopic 
redshifts for 633 galaxies from the Keck Team Redshift 
Survey (KTRS, IWirth et all |2004[ ) were supplemented 
with 695 determin ed using photometric techniques (see 
iBundv et al.ll200^ . The full sample of morphologically- 
selected spheroidals includes 393 sources. Using deep 
Ks-hand imaging obtained at Palomar Observatory for 
GOODS-N and pubhcly available data in GOODS- 
S, stellar masse s were estimated assuming a Chabrier 
()Chabrieil [2001 initial mass f unctio n (IMF) using the 
code described in' Bundv et al.l ((2006')'*. For galaxies with 
photometric redshifts, an additional 0.3 dex is added in 
qua drature to the base uncert ainty which is typically 0.2 

dex (|Bundv et al.l[2005l. [200^ 

As discussed in|Bunds_eLaI| (|2005( ). the GOODS mor- 
phological sample is complete in the JTs-band but be- 
comes incomplete for very red galaxies below a certain 
stellar mass limit as a result of the zab < 22.5 limit (see 
their Figure 4). By considering synthetic SEDs repre- 
senting red (in z — Kg color) stellar populations, we de- 
termine mass limits oi\og{Miimit/MQ) — [9.8, 10.4, 10.8] 
for redshifts, z = [0.4,0.7,1.0]. The completeness lim- 
its of the MFs below are determined at the near edge 
of each redshift interval. When comparing to Mdyn 
we employ the more strict, volume-limited limits at the 
far ed ge of each interval (see discussion in iBundv et al.l 
I2006D . 

3. LINKING DYNAMICAL MASSES AND STELLAR MASSES 

We begin by restricting our analysis to those 125 
spheroidals in GOODS-N for which both Mdyn and M* 
measurements are available. The aim is to develop scal- 
ing relations that can be apphed to the larger GOODS- 
N/S sample. For local galaxies, a fairly tight correlation 
of the form Mdyn = (M*)" has been observed by several 
authors (e.g.,[P admanabhan ct al. 2004; CappcUar i et al.l 
I2006t iGallazzi et al.. ,2006) . .Gallazzi et ah' (,2006[ ) ex- 
ploited a local sample of early-type systems selected by 
concentration, and find a — 1.28 ± 0.03 — the zeropoint 
is not explicitly solved for. In this work, we adopt a 
zeropoint su ch that Mdvn = M<. at lO^^Mo, and also ap- 
ply it to the IGallazzi et al.l (|2006[ ) relation. This choice 
corresponds to Ky = 4.8 and 6.6 at 0.4 < z < 0.7 and 
0.7 < z < 1.0 and is motivated by the fact that the 
depends on the mass density profile and, most impor- 
tantly, on the spectroscopic aperture used to measure 
(J. Since the aperture is effectively redshift-dependent 
(the slit width is fixed at 1"), assuming a constant Ky 
could introduce a spurious redshift dependency. Simi- 
larly, while there is no evidence that t he density profile 
evolv es for high mass spheroidals (e.g.. iKoopmans et al.l 
l2006f) , potential evolution could mimic changes in the MF 
if Ky were held fixed. More detailed studies are needed 
in order to determine the normalization independently. 

We examine the Mdyn-M^, relation at 0.4 < z < 0.7 
and 0.7 < z < 1.0 in Figure [U restricting the compari- 
son to systems above our completeness limit; the best fit 

* The results can be adjusted to those for a Salpeter IMF by 
adding ~0.25 dex to the derived masses. 
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Fig. 1. — Correlation between Mdyn and M«. The dotted line 
illustrates the fitted slope, a, which is also indicated in the top 
lef t corner. A typical error bar is shown and the relation found 
bv lGallazzi et aU | |2006|) is also plotted. The normalization in each 
panel has been set to the same scale, as described in the text. 

values for a are also indicated. O ur observed slopes are 
formally consistent wit h those of IGalla zzi et al.l ()2006l . 
1.28±0.03) at z = and[Rettura et al. (200i, 1.25±0.05) 
at z 1 (|di Serego Ahghieri et al. .2005. , see also), al- 
though some differences may be expected as a result of 
sample selection and fiber- versus slit-based measures 
of velocity dispersion. The fact that a > 1 reflects 
the "tilt" of the Fundamental Plane (e.g. iFaber et al.l 
119871; ICiotti et al.lll996l; iRenzinil 120061) and is consistent 
with non-homology (e.g. iTruiillo et al]|2004D including 
an increasing dark matter f raction with i n Rp at higher 
masses (e.g. [Padm anabha n et al.l 12004 IGallazzi et al.l 
120061 iBolton eral.ll2007f l. With regard to what follows, 
the key result in Figure [1] is that a shows no measur- 
able evolution across our sample. In gl we will show 
that this leads to dynamical MFs that behave similarly 
to their stellar mass counterparts. For consistency, final 
Mdyn estimates for galaxies in the full GOODS sample 
are determined using the value of a determined for the 
two redshift bins, and the final uncertainty in Mdyn is 
taken to be the sum in quadrature of the uncertainty in 
M* and the scatter in Figured] 

4. THE DYNAMICAL MASS FUNCTION 

Using the Mdyn estimates derived above for the full 
GOODS sample, we now construct the dynamical MF 
and compare with its stellar mass counterpart in the top 
two panels of Figure |2l We use the l/V max estimator 
and cl osely follow the method described in IBundv et al.l 
(|2005l ). The solid points with grey shading show the 
stellar MF for the full sample, while the crosses indicate 
the full stellar MF that results when the magnitude limit 
is relaxed to zab < 23.5. The bottom panel of Figure 
121 shows the differential growth derived across the two 
redshift bins. 

Figure 13 demonstrates consistency between the stellar 
MFs and dynamical MFs. Although our adopted zero- 
point in the Mdyn relations ensures a similar vertical nor- 
malization, the application of our scaling relations does 
not mean both MFs should have the same shape. Indeed, 
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Fig. 2. — Comparison of evolving dynamical and stellar mass 
functions. The top two panels show the MFs in two redshift in- 
tervals, while the bottom panel illustrates the differential growth 
between them. Solid points with gray shading trace the stellar 
mass function of the full (all types) GOODS zab < 22.5 sample, 
while crosses denote the equivalent zab < 23.5 sample. Aster- 
isks with light red shading trace the stellar MF of the spheroidal 
component, while diamonds with blue shading indicate the dy- 
namical MF. The fact that the observed increase with time in the 
spheroidal MFs is more significant than the increase observed for 
the full sample (bottom panel), shows that this result is robust 
to potential incompleteness. Isolated error bars indicate the esti- 
mated uncertainty from cosmic variance. 

small difTerences are noticeable in our highest mass bin. 
The good agreement results in large part from the fact 
that the M^yn-M^ relation does not vary significantly 
across our sample. Had a variation in a with redshift 
been present in Figure [U the resulting dynamical MFs 
would evolve with respect to the stellar MFs. The fact 
the two agree validates earlier stellar mass estimates at 
z I and the top-down growth in the abundance of 
spheroidals derived from those measures. 

5. TESTING THE ROLE OF MERGERS 

We now turn to the key question which our analysis can 
address: is the rising abundance with time of spheroidal 
galaxies consistent with that predicted by the growth of 
merging dark matter halos? To answer this question we 
adopt a simple, observationally-motivated model to con- 
nect a and to the virial mass, M^r, of the host ha- 
los. Lensing studi es of early-type galaxies at z w 0.2 
iGavazzi et al.l (see l2007l and references therein) show 
that the total mass density profile of spheroidal halos 
is nearly isothermal within the halo scale radius, r^. We 
therefore assume a mass profile that is isot hermal for r < 
Ts an d follows a p oc r^'^ profile (NFW, iNavarro et al.l 
Il997( ) for r>rs. We note that the p CK r behavior at 
large radii (r ^ r^r) is consistent with non-NFW scal- 
ing laws which can disagree instead at small radii (e.g., 
iGraham et ani2006f ). The mass profile assumed here is 
given by p{r) = a^[2nGr^{l + r/'fTs)]^^, where a is the 
observed central velocity dispersion and we have intro- 
duced the scaling parameter, 7, which can be used to 
tune the profi le shape to observations. By comparing to 
IGavazzi et al.l ( |2007 ), we find 7 — 12 — future observa- 
tions will constrain 7 as a function of redshift and mass. 



We define Myir = ^A{z)pcC'^r^, where pc is the cosmic 
critical density, A(z) i s the overdensity parameter (see 
iBrvan fc Normaiilll998[ ). and r„ir = cr^. The concentra- 
tion, c, is assumed to the follow the relation found in sim- 
ulations, c = 9(1 -f z )-VM„,>/(8.12 x lOi^/i-iMo))-^ " 
(|Bullock et al.ll200ll ). Using this model and the relations 
above, we solve for Myir as a function of velocity disper- 
sion and find, 

l0gM„^ = A^log(T200 +S^+C^(l + ^), (1) 

where M^r has units of HjqMq, (T200 is the velocity dis- 
persion in units of 200 km s~^, and the fit parameters 
are given by Aa ~ 3.1, = 13.2 and Co- = —0.3. Equa- 
tion [1] follows from the model described above and so 
random uncertainties on the fitting parameters are neg- 
ligible. However, systematic errors resulting from tuning 
this mod el to ob servations should refiect the uncertainty 
found bv lGavaz zi et al. (2007) of 0.2 dex. Using our kine- 
matic sample for which we measure both a and M*, we 
can also directly compare M^r and M*, finding the fol- 
lowing relation: 

log My,r = ^* log + + (1 + z) (2) 

where a Chabrier IMF was assumed as before and M^n 
has units of lO^^h^^MQ. The best fitting parameters 
are A = 1-0 ± 0.2 and = 12.83 ± 0.11 ± aB,sys, 
C* = —0.2 ± 0.1 where as.sys = 0.25 dex and ac- 
counts for the systematic uncertainty in Equation [1] (0.2 
dex) as well as uncertainty in the stellar IMF (0.15 
dex). Equations 1 and 2 are applicable for galaxies with 
logM^/Mg « 10 - 12 and cr « 90 - 300 kms'^ Equa- 
tion [2] gives values of Myjr/M ^ = 30l^5, con sistent with 
iMandelbaum et al.l (|2005D and lGavazzi et all (|2007t ). We 
note that our marginally non-zero value of C», if con- 
firmed, implies a slight increase with time of Myir/M^ 
(^30%) over 0.4 < z < 1 .0 that is co nsistent with halo 
growth (see discussion in IConrov et al. 2007) and inde- 
pendent evidence that stel lar mass assembl y is mostly 
completed by z 1 (e.g. iTreu et all |2005( ). We also 
note that virial mass and stellar mass turn out to be ap- 
proximately in linear proportion within this mass range. 
Studies with a larger dynamical range in stellar masses 
and velocity dispersion, and more weak lensing measure- 
ments, are needed to measure the slope with higher pre- 
cision and detect departures from linearity, expected as a 
result of varying star formation efficiency with halo mass. 

Using Equation [2] we estimate halo virial masses for 
our sample and compare the resulting virial MFs in the 
top panel of Figure [31 The dashed line is the total 
halo mass function predicted at z « 0.5 (it evolves little 
over < z < 1. 0) using the online Millennium Simula- 
tion Database^ (jSpringel et a l. 2005b) and adjusted in 
number density by -1-0.1 dex to be consistent with the 
spheroidal MFs (this adjustment is within the expected 
cosmic variance of 0.15 dex). The dotted line shows a lo- 
cal (z w 0.1) estimate constructed by converting the 
MF for early -types (bas e d on concentration, C > 2.6) 
measured by iBell et all (|2003t) into a virial MF using 

^ http://www.mpa-garching.mpg.de/millennium. The "milli- 
Millennium" Database covers a volume of 7 X 10^ h^^Mpc^, nearly 
4 times larger than our largest sampled volume at 0.7 < 2 < 0.4. 
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Equation[2] Adjustments of +0.1 dex were also applied 
to the local number density and the local mass scale to 
match the GOODS sample at high masses. Such cor- 
rections are likely needed as a result of cosmic variance, 
differing selection and mass estimate methods, and pho- 
tometric uncertainties. Thus, comparisons between our 
results and the local MF must be done with caution. 

The top panel of Figure [3] reveals a similar pattern of 
evolution as seen in Figure [5) A useful interpretation 
of this evolution is given in the bottom panel of Figure 
[3] which plots the log difference (log AdN/ dM) , thereby 
displaying the halo MF of new spheroidals that have ap- 
peared over the time interval separating the two redshift 
windows and providing a way to compare with the in- 
cidence of halo mergers. We use the MPAHalo milli- 
Millennium database to identify halos that have expe- 
rienced a recent merger, at which point two previously 
separate halos or subhalos become a single halo. Halo 
merging often occurs several Gyr after halo accretion, 
defined by the point at which the two progenitor halos 
first become associated with the same friends-of-friends 
(FOF) group. Because the smaller progenitor halo may 
experience significant tid al stripping after halo accretion 
but before merging (see iGao et al.l l200l . we take the 
merger mass ratio to be the initial ratio before halo ac- 
cretion. 

After integrating over the two time intervals that sepa- 
rate our observations, we overplot two merging scenarios 
in the bottom panel of Figure[2] (the +0.1 dex vertical off- 
set is also applied). It is a striking result that major halo 
mergers (mass ratios greater than 1:3, indicated by dark 
lines) cannot account for the rise of spheroidal systems 
over the two epochs probed here. The light yellow lines 
(without symbols) show this is still the case when a ~0.5 
Gyr delay is added to the merger timescale (as suggested 
by expected differences in the galaxy and halo merger 
times) . We note that at the highest masses probed there 
is an indication that the predicted merger frequency rises 
above the growth of new spheroidals. This might be ex- 
pected since most systems at such masses already ex- 
hibit spheroidal m orphologies by these redshifts (e.g., 
iBundv et al.ll2005l ). sug gesting such interac tions repre- 
sent "dry" mergers (e.g.. Ivan Dokkuml[200l) . 

6. DISCUSSION 

By constructing dynamical mass functions for 393 field 
spheroidals and linking those mass estimates to the ha- 
los in which they reside, we have found a significant 
result in the bottom panel of Figure [3J major merg- 
ing of dark matter halos as described by ACDM does 
not occur frequently enough to explain the observed in- 
crease in the mass function of spheroidal galaxies. A 
check made by appealing to the semi-analytic model of 
Ide Lucia fc Blaizoti (|2007l ) shows a similar result when a 
galaxy's stellar mass is used to define recent mergers and 



even when is used as the "accounting variable." 

If merging is not the only mechanism that produces 
newly-formed spheroidals since z ^ 1 , other physical pro- 
cesses likely play a significant role. These may include 
mechanisms for transforming disk and irregular galax- 
ies into relaxed spheroidals, for example secular bulge 
growth accompanied b y disk fading. Indeed , recent semi- 
analytic models fe.g.. [Croton et all 120061 : iBower et al.l 
120061: Ide Lucia fc BlaizotI 120071) have required prescrip- 
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Fig. 3. — Evolving virial mass functions of halos hosting 
spheroidals. The top panel shows MFs in two z-bins. The local 
MF estimated fromiBell et al. (2003), is indicated by the dotted 
line. The dashed line shows the total MF of predicted dark matter 
halos, and the isolated error bar (top right) denotes the cosmic 
variance uncertainty. The bottom panel shows the MF of newly 
formed spheroidals determined by subtracting the MFs (top panel) 
in adjacent redshift bins. Overplotted are two halo merger scenar- 
ios based on the Millennium Simulation, integrated over the red- 
shift ranges indicated. Black lines with triangle symbols trace >1:3 
mergers while yellow lines (without symbols) denote such mergers 
with a ~0.5 Gyr time delay. 

tions for such mechanisms, in addition to mergers, to 
match the local abundance of spheroidals. More work is 
needed to understand these processes and their role in 
morphological evolution. 
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